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INTRODUCTION 


This  report  describes  a  solution  foraulatlon  for  and  Its  applications  to 
Initial  boundary  value  probleas  of  structural  dynamics  and  stress  waves. 
Excellent  numerical  results  are  stated  In  conjunction  with  finite  eleaent 
discretization.  The  basic  concept  of  this  approach  Is  to  establish  a  varia¬ 
tional  problea  equivalent  to  a  given  Initial  boundary  value  problem ,  which  Is 
In  general,  non-self-adjolnt,  through  the  use  of  an  adjoint  field  variable  and 
the  use  of  soae  large  "spring"  constants  so  that  all  the  end  conditions  can  be 
transformed  into  natural  "boundary"  conditions.  Therefore,  the  shape  func¬ 
tions  used  need  not  satisfy  any  end  conditions  a  priori  In  solving  the  varia¬ 
tional  problea  In -the  same  manner  as  applying  the  Raylelgh-Ritz  method  for 
self-adjoint  probleas.  This  saae  concept  was  demonstrated  In  solving  Initial 
value  probleas  in  a  paper  delivered  at  the  International  Symposium  on 
Numerical  Methods  In  Engineering  Science  series  in  1978  and  later  published  in 
the  Journal  of  Sound  and  Vibration. 1  In  this  present  report,  the  foraulatlon 
is  extended  to  Initial  boundary  value  probleas  and  the  numerical  results 
obtained  are  also  encouraging. 

In  the  section  which  follows  immediately,  two  Initial  boundary  value 
probleas  are  stated.  One  is  a  longitudinal  stress  wave  problea  In  a  rod. 
There  Is  a  discontinuity  In  the  Initial  data  given.  He  wish  to  trace  this 
discontinuity  In  the  numerical  solution  using  the  present  apprr«ch.  The 
second  problea  Is  a  beam  vibration  problea  under  a  moving  concentrated  load. 

lj.  J.  Hu,  "Solutions  to  Initial  Value  Probleas  by  Use  of  Finite  Elements  - 
Unconstrained  Variational  Formulations,"  1977  Journal  of  Sound  and  Vibration, 
33,  pp.  341-356. 


This  is  a  auch  more  difficult  problem  alnca  the  partial  differential  equation 
la  of  fourth  order  and  the  force  la  singular  In  nature.  In  the  next  section, 
variational  probleas  equivalent  to  the  given  Initial  boundary  problems  are 
established.  The  finite  eleaent  discretisation  procedures  are  then  briefly 
recaptured.  Lastly,  nuaerlcal  results  are  preaented  with  soae  comments. 

INITIAL  BOUNDARY  VALUE  PROBLEMS 

Two  Initial  boundary  probleas  of  structural  dynamics  will  be  stated  In 
this  Section.  The  first  one  is  of  longitudinal  elastic  stress  wave  In  a  rod 
with  a  sudden  change  in  Initial  conditions.  The  second  one  Is  concerned  with 
lateral  vibrations  of  a  Euler -Bernoulli  beam  subjected  to  a  moving 
concentrated  load. 

Longitudinal  Stress  Wave  in  a  Rod 

The  rod  Is  fixed  at  one  end  and  free  at  the  other  end.  The  diacontlnulty 
data  arises  from  the  Initial  linear  displacement,  corresponding  to  a  constant 
stress,  due  to  a  force  applied  at  the  "free"  end.  This  force  suddenly 
disappears  at  tlaa  sero  causing  a  stress  discontinuity  at  the  free  end.  The 
differential  equation  can  be  written  as: 

3*u  1  32u  0  <  x  <  1 

rr  "  *7  T?  *  U> 

3m*  a 1  3tz  0  <  t  <  T 

with 

a2  -  E/p  (2) 

where  u  ■  u(x,t)  is  the  axial  displacement ;  x,t  are  the  coordinates  In  axial 
direction  and  In  tlaa,  respectively;  p,B  are  density  and  Young's  modulus, 
respectively ,  of  the  rod  material ;  1  denotes  length  of  the  rod ;  and  T  denotes 
soaa  finite  tfaw  of  Interest. 
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For  Che  boundary  condition*,  we  hare 


and 


u(0,t)  »  0 

3u  <3) 

~  (I.t)  -  0 


Th*  dynaalcs  of  the  problea  are  due  to  the  initial  conditions.  It  is  assuaad 
that  the  rod  is  stretched  to  a  linear  dlsplaceaent  by  a  force  P  which  vanish** 
at  tlae  t  >  0  (see  Figure  1).  The  initial  velocity  of  the  rod  is  assuaad  to 


Figure  1.  A  Bod  Fixed  at  One  Bad  and  Subjected  to  a  Load  P,  which  is 
Suddenly  Beleased  at  Tlae  Zero. 

It  is  convenient  to  use  dlaenslonlees  parsaeters.  Let 

u*  •  u/i  ,  &  “  x/1  ,  t*  •  t/T  (5) 

Then,  Bq.  (1)  In  dlaenslonlese  fora  is 

32(1*  „  0  <  x*  <  1 

“TJ  -  b2  —r  ,  (6) 

*x*2  3t*2  0  <  t*  <  1 
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vh«r« 


Tht  boundary  conditions  become 


o  1  1 

b  -  1  <:> 
a*  T 


au* 


and 


u*(0,t*)  -  0  ,  (l,t*)  -  0 

«• 


du* 

u*(x,0)  “  P*x*  ;  -  (x*,0)  •  0 

at* 


where 

P 

P*  •  — 
AS 


(7) 

(8) 

(9) 

(10) 


is  the  force  In  dimensionless  fora* 

The  stated  problem  in  dlaenslonless  fora  coablnes  Eqs.  (6) ,  (8) ,  and  (9) 
with  the  new  dlaenslonless  parameters  related  to  physical  counterparts  by  Kqs* 
(5),  (6),  and  (10).  To  slapllfy  writing,  as  shall  drop  the  asterisks  (*)  In 
Bqs.  (6) ,  (8) ,  and  (9) ,  and  rewrite  thaa  as 


- 

0  <  x  <  1 

u"  -  b2u  ■  0 

• 

• 

(6*) 

0  <  t  <  1 

u(0,t)  -  0  ; 

u'(l,t)  -  0 

(8*) 

u(x,t)  *  Px  ; 

u(x,0)  "  0 

(9‘) 

where  a  prlaa  (')  indicates  differentiation  with  respect  to  x  and  a  dot  (•) , 
with  respect  to  t. 

Baaa  Vibrations  Under  Mowing  Loads 

Let  us  consider  the  differential  equation  of  a  uniform  Euler-Bernoulll 
beam  subjected  to  a  aorlng ,  concentrated  force. 

0  <  x  <  1 

Ur""  ♦  pAy  -  P«(xp-x)  (11) 

0  <  t  <  T 


4 


where 

E,p  “  Young's  aodulus,  density  of  the  besa  asterlal 

I,A  ■  second  aoaent,  ares  of  the  bees' s  cross-section 

1  ■  length  of  the  bean 

y*y(x,t)  ■  beam  deflection 

x,t  »  coordinates  In  besa's  axial  direction  and  In  tiae 

P  ■  aagnltude  of  the  concentrated  force 

((x)  •  Dirac  delta  function 

Xp-Xp(t)  “  location  of  P 

T  -  soaa  finite  tlae  of  Interest 
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VARIATIONAL  PROBLEMS  -  A  GENERALIZED  RAYLEIGH-RITZ  METHOD 

For  the  stress  wave  problem  In  the  previous  section,  consider  s 
verlstlonsl  problem. 

6lo  -  0  (16s) 

with 

,1.1  .  • • 

I0  •  Io(u,v)  -  J  (-u'v,+b2uv)d*dt  (16b) 

where  u(x,t)  and  v(x,t)  are  said  to  be  adjoint  to  each  other.  It  Is  a  simple 
matter  to  see  that  this  problem  Is  an  Indeterminate  one.  However,  the 
functional  of  Eq.  (16b)  can  be  modified  to  a  variational  problem  which  Is 
equivalent  to  the  boundary/initial  problem  of  Eqs.  (6'),  (8'),  and  (9').  Thus 
consider 

61-0  (17a) 

with  11 

I  -  I(u,v)  -  /  /  (-u'v'+b2uv)dxdt 

0  0 

+  M/*  u(0,t)v(0,t)dt 
0 

„,1  ,  ,1 

♦  *2b 2/q  [u(x,0)  -  Uo(x)]v(x,l)dx  +  b2  /  ui(x)v(x,0)dx  (17b) 

He  shall  take  the  first  variation  of  the  function  I(u,v)  of  Eq.  (17b)  In 
such  a  manner  that  6v  Is  completely  arbitrary  while  6u  Is  set  to  sero 
Identically.  Hence,  by  means  of  lnt egr at lon-by-par t s ,  one  has 
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(6I)6u»0  "  /  /  ( u“ -b2u) Svdxdt 

0  0 

-  J  u'(l,t)6v(l,t)dt 
0 

» 

+  J  (u(0,t)  +  k ju(0 ,t) ] dv( 0 ,t  )dt 

+  k2/  (u(x,l)  ♦  k2(u(x,0)  -  Uo(x)]}(v(x,l)dx 

0 

,  ,1  . 

-  b2  J  (u(x,0)  -  ux(x)]6v(x,0)dx  »  0 
0 


The  face  that  6v(x,t)  la  completely  arbitrary  aaaklaa  ua  to  conclude  frou  Bq. 


(18)  that 


u"  -  h2u  »  0  ; 


0  <  x  <  1 


0  <  t  <  1 


(19a) 


u'(l.t)  -  0 
u'(0,t)  +  kxu(Ott)  •  0 
u(x,l)  ♦  k2lu(x,0)  -  Uo(x)]  -  0 
u(x,0)  -  ux(x)  -  0 


(19b) 


It  la  than  obaarrad  that  tha  lnltlql  boundary  value  problaa  daflnad  by  Iqa. 
(19a)  and  (19b)  raducaa  to  that  of  Iqa.  (6'),  (8*).  and  (9')  If  ona  lata  kx 
and  k2  go  to  Infinity*  (and  with  Uq(x)  -  Px  and  ux(x)  •  0).  Thla  fact 
auggaata  that  tha  variational  problaa  of  Iqa.  (17a)  and  (17b)  can  ba  uaad  aa  a 
baala  of  a  flnlta  alaaant  dlacratlaatlon  for  tha  approxlaata  aolutlona  to  tha 
original  Initial  boundary  problaa.  It  ahould  ba  notad  that  all  tha  auxiliary 


*1hla  procaaa  la  aoaatlaao  rafarrad  to  aa  tha  paaalty  function  aathod.  See, 
for  axaapla,  lafaranea  2. 

*D.  0.  Luanbergar,  Oatlalaatlon  by  factor  Spaca  Hithod.  John  HLley,  1969,  p. 
302. 


conditions  In  Eqs.  (19a)  and  (19b)  are  the  so-called  natural  boundary 
conditions.  They  are  the  consequence  of  the  variational  problem  -  just  like 
the  differential  equation  itself,  for  this  reason,  the  above  solution  Is 
referred  to  as  a  Generalised  Raylelgh-Rlts  Method. 

By  a  similar  process,  one  can  establish  a  variational  problem  for  the 
vibration  problem  of  a  beam  under  a  moving  load.  In  this  case,  one  has 

.1  ,1  .  .  -  - 

61  ■  /  /  [u"6v"  -  u6v  -  6(x-x)6v]dxdt 

0  0 

rl 

+  /  lkiu(0,t)6v(0,t)  +  k2u’(0,t)6v'(0,t) 

0 

.  +  k3u(l,t)6v(l,t)  +  k4u'( 1 ,t)6v'( 1 ,t) Jdt 
,1 

+  /  [k5u(x,0)6v(x,l)  +  k£u(x,0)6v(x,0) ]dx  -  0  (20) 

Through  lntegratlons-by-parts, 

11  ~ 

61  ■  /  /  lu""  +  u  -  6(x-x) ]6v(x,t)dxdt 

0  0 

+  I1  {[k1U(0,t)  +  u“'(0,t)]6v(0,t)  +  [k2u'(0,t)  -  u“(0,t)]6v'(0,t) 

+  (k3u(l,t)  -  u-»(l,t)J6v(l,t)  +  [k4u'(l,t)  +  u-(l,t)]«v*(l,t)}dt 

+  /  {(k5(u(x,0)-0)  -  u(x,t)]6v(x,l)  +  (k^+l) Iu(x,0)-0l 6v(x,0) }dx  -  0  (21) 

0 

The  original  differential  equation  and  the  boundary  and  Initial  conditions  are 
recovered  from  the  equation  above  due  to  the  arbltrarlnees  of  the  variations 
6(x,t)  and  by  properly  selecting  the  values  of  kj.,s,  1  *  1,2,.. .,6. 
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FINITE  ELEMENT  DISCRETIZATION 


Only  essential  feature*  will  be  stated  In  the  finite  eleaant  discretisa¬ 
tions  here.  The  region  of  a  unit  square  (0  <  i  <  1;  0  <  t  <  1)  Is  further 
divided  into  KxL  rectangles  by  taking  K  divisions  In  x  direction  end  L 
divisions  In  t  direction.  Local  coordinates  (£,n)  In  each  (i,J)th  eleaant  are 
related  to  (x ,t)  by  these  equations* 

t  -  e<4>  -  to  -  1  ♦  l 

.  (22) 
n  -  n(^  ■  Lt  -  j  +  l 

Within  each  eleewnt ,  the  unknown  function  u(x,t)  Is  replaced  by  the 
approximation: 

u(l,J)(5.n)  -  sTCC.»i)  0(t. j) 

(23) 

*r(i,j)<*»l'>  "  •*<*»*>  *v(i,j) 

where  a(€,ri)  is  the  shape  function  vector  end  O(i.j).  *V(l,j)  *H  tOe 
generalised  coordinates.  The  specific  fora  of  a(€,n)  saployed  here  Is  such 
that  each  one  of  the  sixteen  coaponents  1st 


•k(C.n)  -  bi(Obj(n)  . 


k  -  1.2, 


1.3  -  1. 2.3.4 


bi(e)  -  1  -  3C2  +  2C3 
b2(0  -  K  -  2tz  +  C3 

I 

h3(»  -  3C*  -  2C3 
b4(0  -  -C2  +  5  3 

and  the  relations  between  Index  k  and  the  pair  (l.j)  are  given  In  Table  I. 


TABLE  I.  RELATIONSHIP  BETWEEN  (l,j)  AND  k  IN  EQUATION  (24) 


k 

a.i> 

k 

U.l> 

1 

(i,D 

9 

(1,3) 

2 

(2,1) 

10 

(2,3) 

3 

(1,2) 

11 

(1,4) 

4 

(2,2) 

12 

(2,4) 

5 

(3,1) 

13 

(3,3) 

6 

(4,1) 

14 

(4,3) 

7 

(3,2) 

15 

(3,4) 

8 

(4*2— 

16 

Using  Eqs.  (22)  through  (25)  In  Bq.  (17)  «nd  ths  fact  that  V(i,j)  Is 
completely  arbitrary,  the  Matrix  aquations  for  tha  unknowns  U(l,j)  can  ba 
routinely  assembled  and  solved,  further  details  will  be  oulttad  here. 

NUMERICAL  RESULTS  AMD  DISCUSSION 

Soae  of  the  numerical  results  -ere  presented  In  this  section.  Ibr  the 
stress  neve  problem*,  table  II  provides  solutions  of  v(x,t),  3u/dx(x,t)  and 
3u/3t(x,t)  for  x  ■  0,  0.1,  0.2,  ...1.0  and  for  t  »  0,  0.5,  1.0,  1.5,  and  2.0. 
During  this  tlee  Interval,  the  original  dlsplacemant  has  gona  through  a 
complete  sign  reversal  as  shown  In  Figure  2.  This  particular  set  of  data  was 
obtained  by  taking  K  ■  10  and  L  •  1  with  restart  procedures,  l.e.,  the  final 

*#or  exact  solution  to  this  problem,  see  for  example.  Reference  3. 

^L.  S.  Jacobson  and  R.  S.'  Byre,  Engineering  Vibrations,  McGraw-Hill ,  1965,  pp. 
472-474. 
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solution  in  the  first  time  step  was  taken  as  the  initial  condition  of  the  next 
step  in  time,  and  so  on*  Values  of  the  exact  solutions  are  given  in 
parentheses*  Excellent  agreement  is  observed*  The  fact  that  the 
discontinuity  of  the  solution  follows  along  without  much  oscillation  is  worth 
mentioning* 

For  the  beam  vibration  problem  with  a  moving  force,  some  typical 
numerical  solutions  are  given  in  Tables  III  and  IV.  The  moving  concentrated 
force  is  assumed  to  travel  at  a  constant  velocity  c  (although  this  Is  not  at 
all  a  restriction  for  the  present  method)  such  that 

x(t)  -  ct 

where  c  is  dimensionless  velocity.  For  small  c,  c  -  0.0001,  and  the 
displacement  solutions  become  those  of  static  deflections  as  shown  in  Table 
III.  For  a  large  c  (compared  with  unity),  c  *  10,  and  solutions  show  dynamic 
effects  as  indicated  in  Table  IV*  As  a  comparison,  solutions  obtained  by  the 
Fourier  series  and  Laplace  transform  method^  are  given  in  parentheses*  Good 
agreement  exists  even  in  cases  with  considerable  dynamic  effect. 

In  conclusion „  this  report  has  demonstrated  through  examples  of 
structural  dynamics  an  approximate  solution  formulation  (which  is  both  a 
weighted  method  and  a  variational  problem),  the  finite  element  implementation, 
and  some  favorable  numerical  results*  Although  only  linear  problems  have  been 
mentioned,  an  application  to  solutions  of  non-linear  problems  is  now  being 
investigated. 

*L.  Fryba,  Vibrations  of  Solids  and  Structures  Under  Moving  Load,  Noordhoff, 
1971. 


TUU  III.  SOLUTIONS  u(*,t)  TO  THE  NOTING  PORCH  PROBLEM  OF  BQ.  (15) 


TAILS  IV.  SOLUTIONS  u(x,t)  TO  THE  MOVING  FORCE  PROBLEM  OF  EQ.  (15) 
WITH  Q  -  1.0  AND  FIXED  END  CONDITIONS  AT  x  -  0. 

(For  very  low  velocity,  Y  ■  10  in  Eq,  (14)) 
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